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We report results of theoretical studies on the elastic properties of single-wall nanotubes of the 
following compositions: C, BN, BC3, BC2N and C3N4. These studies have been carried out using a 
total energy, non-orthogonal tight-binding parametrisation which is shown to provide results in good 
agreement both with calculations using higher levels of theory and the available experimental data. 
Our results predict that of all types of nanotubes considered, carbon nanotubes have the highest 
Young's modulus. We have considered tubes of different diameters, ranging from 0.5 to 2 nm, and 
find that in the limit of large diameters the mechanical properties of nanotubes approach those of 
the corresponding flat graphene-like sheets. 



The discovery of Ceo and fullerenes 0] in the mid 80 's 
was soon followed by the observation of nanotubes 0, 
first reported by Iijima || in 1991. Since then nanotubes 
have been the focus of attention of a growing scientific 
community, attracted to them by their many interesting 
properties, such as their structure, electrical conductiv- 
ity and mechanical properties, as well as by their large 
potential for practical applications. Two types of nan- 
otubes exist: those originally observed by Iijima || were 
multi-wall nanotubes (MWNT's), formed by concentric 
shells of apparently seamless cylinders of graphcnc, hav- 
ing a separation between them similar to that in graphite. 
More recently, single- wall nanotubes (SWNT's) have also 
been synthesized. As their name indicates, these consist 
of a single seamless cylinder of graphcnc . 

Soon after the discovery of carbon nanotubes it was 
proposed that other compounds forming graphite-like 
structures, such as BN {§, BC 3 §, BC 2 N §, and 
CN [7], could also form nanotubular structures. Indeed 
BN 5-0, BC 3 and BC 2 N have now been synthe- 
sized, though the actual structure of BC 2 N tubes seems 
to correspond to concentric shells of C and BN in a 'sand- 
wich' structure |0|. Other tubular structures formed by 
heavier element compounds have been predicted, such as 
GaSe @, and synthesized, like WS 2 and MoS 2 @. 

In this paper we focus our attention on the structural, 
energetic and mechanical properties of single-wall car- 
bon and composite nanotubes. We perform a systematic 
study of these systems using Tight-Binding total energy 
methods |l5|], as well as first-principles Jll| calculations. 
Some of the results reported here have already appeared 
in published form p7| , but we also provide previously 
unpublished results for the C3N4 nanotubes. 

The structure of this paper is as follows. In Sec. | wc 
review the available experimental data on the mechanical 
properties of nanotubes, while Sec. || is devoted to dis- 
cussing previous theoretical work. In Sec. [II we describe 
the models used in our work and the calculations carried 
out in order to address the issues discussed here. Then, 



in Sec. IV we discuss our results and conclusions. 



I. REVIEW OF EXPERIMENTAL RESULTS 

There is a growing body of experimental evidence 
indicating that carbon nanotubes (both MWNT's and 
SWNT's) have extraordinary mechanical properties. 
There are many direct observations of the large bending 
flexibility |^8|-po| of nanotubes, which provide evidence of 
their capability to sustain large strains without evidence 
of collapse or failure. However, the technical difficulties 
involved in the manipulation of these nano-scale struc- 
tures makes the direct determination of their mechanical 
properties a rather challenging task. In spite of these dif- 
ficulties, a number of experimental measurements of the 
Young's modulus of nanotubes have been reported. The 
first such study was that of Treacy et al. |pl[ , who corre- 
lated the amplitude of the thermal vibrations of the free 
ends of anchored nanotubes as a function of temperature 
with the Young's modulus. Regarding a MWNT as a hol- 
low cylinder with a given wall thickness, one can obtain 
a relation between the amplitude of the tip oscillations in 
the limit of small deflections, and the Young's modulus. 
Having quantified the amplitude of those oscillations by 
means of careful TEM observations of a number of nan- 
otubes, Treacy et al. were able to obtain an average value 
of 1.8 TPa for the Young's modulus, though there was 
significant scatter in the data (from 0.4 to 4.15 TPa for 
individual tubes). Thus this number is subject to large 
error bars, but it is nevertheless indicative of the excep- 
tional axial stiffness of these materials. 

More recently Krishnan et al. 0| have reported stud- 
ies on SWNT's using the same technique. A larger sam- 
ple of nanotubes was used, and a somewhat smaller av- 
erage value was obtained, Y — 1.25 TPa, closer to the 
expected value for graphite along the basal plane. 

This technique has also been used by Chopra and 
Zettl |23| to estimate Y for BN nanotubes. Their results 
indicate that these composite tubes are also exception- 
ally stiff, having a value of Y around 1.22 TPa, very close 
to the value obtained for carbon nanotubes. 

Another way to probe the mechanical properties of 
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nanotubes has been described by Wong et al. [Q], who 
have used the tip of an Atomic Force Microscope (AFM) 
to bend anchored MWNT's while simultaneously record- 
ing the force exerted by the tube as a function of the 
displacement from its equilibrium position, information 
from which the Young's modulus of the nanotube can be 
extracted. Wong et al. have reported a mean value of 
1.28 TPa, which is in good agreement with the previous 
experimental results. Also Salvetat and coworkers ]25| ] 
have used a similar idea, which consists of depositing 
MWNT's on an ultra-filtration membrane. Many tubes 
are then found to lie across the holes present in the mem- 
brane, with a fraction of their length suspended. The tip 
of an AFM is then used to exert a load on the suspended 
length of the nanotube, measuring at the same time the 
nanotube deflection. The mean value of the Young's 
modulus obtained by Salvetat et al. was 0.81 TPa. A 
similar procedure has also been used by Muster et al. [^6) , 
who used an AFM to record the profile of a MWNT ly- 
ing across an electrode array. By assuming a simple Van 
der Waals interaction law between the tube and the sub- 
strate, and regarding the nanotube as an elastic beam, 
the measured profile was found to be consistent with a 
Young's modulus of approximately 1 TPa. 

Other experiments, which as yet have not aimed at the 
mechanical characterization of nanotubes, nevertheless 
hint at other possible ways in which this characterization 
could be carried out. Among these we could cite the 
embedding of nanotubes in resins |27j] , or measuring the 
bending of anchored nanotubes in controlled magnetic 
fields |f. 

All these experiments have contributed to confirming 
that nanotubes, both SW and MW, have indeed excep- 
tional mechanical properties, but there are still ques- 
tions that remain unresolved. The work of Chopra and 
Zettl ||^] indicates that BN nanotubes are close in stiff- 
ness to carbon nanotubes, but the experimental error 
bars are too large to state categorically that one type 
of tube is stiffer than the other. The results of Salve- 
tat et al. (||] clearly indicate that MWNT's synthesized 
by the arc-discharge method j2{| are much stiffer than 
those produced by the catalytic decomposition of hydro- 
carbons fsofl , which have a Young's modulus in the range 
of 10-50 TPa. This large difference is presumably a reflec- 
tion of the influence of the high density of defects present 
in the structure of the latter tubes, but a detailed quan- 
tification of the influence of defects on the mechanical 
properties of nanotubes is as yet missing. 



II. PREVIOUS THEORETICAL WORK 

The mechanical properties of nanotubes have been ad- 
dressed also by means of theoretical calculations in a 
number of publications [^8 31-37]]. Most of these stud- 
ies have been carried out using empirical potentials, al- 



though tight-binding based models have also been oc- 
casionally used |33|. Though well-tested empirical po- 
tential models exist for carbon-based systems [[38"|-fl0"| , to 
our knowledge no such model exists for the composite 
systems, and therefore these materials have been mostly 
studied using first-principles methodologies QJ^]. How- 
ever, these latter studies have concentrated largely on 
electronic, structural and vibrational properties of the 
composite systems, without addressing the issue of their 
mechanical properties. 

Lu j37j] has reported an extensive study of the Young's 
modulus, Poisson ratio and elastic constants of car- 
bon nanotubes (both SW and MW, as well as ropes of 
SWNT's) using an empirical pair potential. Yakobson et 
al. |m| and Nardelli et al. J35|] have studied the behaviour 
of nanotubes subject to large axial strains. 



III. MODEL AND CALCULATIONS 

For the majority of the calculations reported here we 
have used a non-orthogonal Tight-Binding scheme due to 
Porezag and coworkers 0]. Tight-Binding (TB) meth- 
ods jl5| lie in the centre region of the spectrum of simu- 
lation methods in both computational cost and reliabil- 
ity. Empirical potentials |38| ] are much cheaper to use, 
but their accuracy and reliability is often questionable. 
First-principles methods |l6| ] on the other hand are more 
reliable, but their computational demands are orders of 
magnitude larger than for TB calculations, and often this 
makes their use impractical. 

The TB model of Porezag et al. |4l]] contains two con- 
tributions to the total energy: a so-called band- structure 
energy term, and a repulsive pair-potential. The band- 
structure energy is calculated as the sum of the eigenval- 
ues of the occupied states of a TB Hamiltonian. A matrix 
representation of this Hamiltonian is constructed using a 
minimal basis set consisting of a single atomic-like or- 
bital per atomic valence state. The matrix elements are 
evaluated in the framework of Density Functional The- 
ory (DFT) (normally in the Local Density Approxima- 
tion, LDA), but retaining only two-centre contributions 
to the integrals. This means that each matrix element 
of the Hamiltonian depends only on the relative distance 
of the two atoms on which the corresponding basis func- 
tions are centred and the direction cosines of the internu- 
clear vector |42|Jl5[| . Since the basis set is not orthogonal, 
the calculation of the band structure energy requires the 
solution of a generalized eigenvalue problem [^3). The 
repulsive pair potential is then constructed in such a way 
that the total TB energy of a reference system (usually 
the dimer) matches that of the full DFT calculation with 
the same basis set. For more details on the construction 
of the model the reader should consult the original ref- 
erences 1 41 , but an important point worth emphasizing 
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here is the fact that no information concerning the me- 
chanical properties of the system under study are used 
in the construction of the TB parametrisation. 

Even thought the TB model used here is not fitted to 
any empirical data, it is nevertheless approximate and 
less accurate than conventional first-principles methods. 
For this reason we have also carried out plane- wave (PW) 
pseudopotential DFT calculations for the (6,6) C and BN 
nanotubes, in order to have the possibility of compar- 
ing our TB results with fully ab initio calculations. The 
PW calculations were performed using Troullier-Martins 
pseudopotentials Q with a PW cutoff of 40 Ry for the 
basis set, and 10 reciprocal space points generated ac- 
cording to the Monkhorst-Pack scheme j|5j to sample 
the one-dimensional Brillouin zone. The TB calculations 
used T-point sampling only, but the periodic cells were 
chosen large enough as to ensure the same degree of con- 
vergence in total energy differences as were achieved in 
the PW calculations. 

As we have seen in Sec. |, the central property charac- 
terizing the stiffness of nanotubes, to which the experi- 
ments have access, is the Young's modulus. In bulk 3-D 
systems Y is given by the following expression: 
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where E is the total energy, e is the strain and V Q is 
the equilibrium volume. The second derivative measures 
how rapidly the energy grows as the system is distorted 
out of its equilibrium configuration. Usually Y is given 
in units of pressure, and this is why the factor of Vq -1 
appears in this formula. However, Eq. (|lj) presents an 
ambiguity in the case of SWNT's, which stems from the 
definition of Vq. To define Vq for a SWNT one needs to 
specify the wall thickness, and there is no clear way to 
define the thickness of a wall one-atom thick. In previous 
publications different authors have used different values 
for the tube wall thickness, thought the most common 
convention has been to adopt the value of the interlayer 
spacing in graphite. Nevertheless the value of Y depends 
on the inverse of the wall thickness 5 R, and is therefore 
rather sensitive to the chosen value. In a recent publica- 
tion fll7|| we proposed a way to bypass this problem by 
using an alternative definition of the Young's modulus, 
more appropriate to the case of SWNT's: 



Y s = 
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Here So is the surface area defined by the nanotube at 
zero strain, which is a well defined quantity. Given that 
Uo = SqSR, one can recover the usual definition by simply 
dividing by 8R: Y = Y S /SR, if one wishes to adopt a 
particular convention. 

Another mechanical property of interest is the Poisson 
ratio, a, which is defined by 
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where R is the radius of the tube at strain e, and R eq is 
the equilibrium (zero strain) tube radius. The Poisson 
ration measures how much the tube contracts (expands) 
radially when subject to a positive (negative) axial strain 
e. 

We have performed a series of calculations using 
the Tight Binding model discussed above aimed at de- 
termining Y s and a of both carbon and composite 
SWNT's,including BN, BC 3 BC 2 N and C 3 N 4 nanotubes. 
Two different graphite-like BC2N structures are possible, 
but in our studies we have only considered the structure 
known as II, since this is the one reported to be most 
stable g§. 

The calculations consist of taking a section of a nan- 
otube using periodic boundary conditions to simulate an 
infinite tube, and subject it to both negative and pos- 
itive strains along the axial direction. At each strain 
the positions of all atoms in the repeat-cell are fully re- 
laxed without constraints, using the Conjugate Gradients 
minimisation technique f43f| . A given calculation was as- 
sumed to have converged once the total energy varied 
less than 10~ 5 Hartree between two successive iterations. 
From these calculations we obtained the total energy and 
atomic positions as a function of the axial strain imposed 
on the nanotube, from which we could calculate Young's 
modulus and the Poisson ratio, as well as the equilibrium 
structures of each nanotube considered. 



IV. RESULTS AND DISCUSSION 

Let us first consider the difference of energy between 
a nanotube structure and the corresponding infinite flat 
graphene sheet. This energy difference is known as the 
strain energy E s , and we have plotted it in Fig. |l] for 
C, BN and BC3 (n,n) nanotubes, as a function of the 
tube diameter. It can be seen from Fig. |l| that the strain 
energy varies as D~ 2 , where D is the tube diameter; we 
have performed fits of functions of the form aD~ b to the 
data, and the values for the parameters a and b are given 
in Table | 
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FIG. 1. Curvature strain energy as a function of the equi- 
librium tube diameter, as obtained from the tight-binding cal- 
culations, for C, BN and BC3 nanotubes. 

That the strain energy should decay as D~ 2 was 
predicted by Tibbetts (see also Mintmire and 

White (48) ) on the basis of continuum elasticity theory, 
according to which the strain energy per atom is given 
by 

E s _Ya 3 fl 

where Y is the Young's modulus of the tube, a is a con- 
stant of the order of the inter-layer spacing in graphite 
and fl is the area per atom. Note that from Eq. (0) and 
the numerical fits to the strain energy data of Fig. [j]given 
in Table |, one can predict that the Young's modulus of 
SWNT's of BN and BC3 of a given diameter should be 
approximately 0.68 and 0.71 respectively that of a carbon 
nanotube of the same diameter. As we shall see below, 
direct calculations of Young's modulus for these tubes 
obey approximately this relation. 



BxCyN z (n,m) a x 10 2 (eV nm 2 /atom) b 

C (n^n) O 2.083 

(n,0) 8.7 1.996 

BN (n,n) 5.5 1.984 

(n,0) 5.6 1.980 

BC 3 (n,n) 5.8 1.984 

(n,0) 5.6 2.048 



TABLE I. Parameters obtained from fitting the strain en- 
ergy curves of Fig. |l| to a function of the form aD~ b . Note 
that the value of b is very close to 2 in all cases. 



B^CyNz 


(n,m) 


D eq (nm) 


a 


Y s (TPa • nm) 


Y (TPa) 


c 


(10,0) 


0.791 


0.275 


0.416 


1.22 




(6,6) 


0.820 


0.247 


0.415 


1.22 






(0.817) 




(0.371) 


(1.09) 




(10,5) 


1.034 


0.265 


0.426 


1.25 




(10,7) 


1.165 


0.266 


0.422 


1.24 




(10,10) 


1.360 


0.256 


0.423 


1.24 




(20,0) 


1.571 


0.270 


0.430 


1.26 




(15,15) 


2.034 


0.256 


0.425 


1.25 


BN 


(10,0) 


0.811 


0.232 


0.284 


0.837 




(6,6) 


0.838 


0.268 


0.296 


0.870 






(0.823) 




(0.267) 


(0.784) 




(15,0) 


1.206 


0.246 


0.298 


0.876 




(10,10) 


1.390 


0.263 


0.306 


0.901 




(20,0) 


1.604 


0.254 


0.301 


0.884 




(15,15) 


2.081 


0.263 


0.310 


0.912 


BC 3 


(5,0) 


0.818 


0.301 


0.308 


0.906 




(3,3) 


0.850 


0.289 


0.311 


0.914 




(10,0) 


1.630 


0.282 


0.313 


0.922 




(6,6) 


1.694 


0.279 


0.315 


0.925 


BC 2 N II 


(7,0) 


1.111 


0.289 


0.336 


0.988 




(5,5) 


1.370 


0.287 


0.343 


1.008 


C3N4 


(6,0) 


0.913 


0.280 


0.192 


0.565 




(8,0) 


1.210 


0.238 


0.207 


0.610 




(6,6) 


1.558 


0.177 


0.228 


0.670 




(8,8) 


2.075 


0.132 


0.233 


0.684 



TABLE II. Structural and elastic properties of selected 
nanotubes obtained from the tight-binding calculations re- 
ported here. Young's modulus values given in parenthesis 
were obtained from first-principles calculations. Also the 
value of Y with the convention SR = 0.34 nm is given for 
comparison. 
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A first indication that the TB model used in this work 
is a reliable one stems from the good agreement obtained 
in the strain energy as calculated here and that calcu- 
lated from first-principles methods and reported else- 
where . Another indication comes from the fact that 
a certain buckling on the surface of the BN nanotubes 
is predicted to occur, also in agreement with preliminary 
first-principles calculations [0]. This buckling, which re- 
sults from the B atoms displacing inwards towards the 
tube axis, while the N atoms displace in the opposite 
direction, is a consequence of the slightly different hy- 
bridizations of B and N on the curved surface of the nan- 
otube. The amount of buckling is dependent on the tube 
diameter, but it is otherwise independent of the tube 
structure for arm-chair and zig-zag nanotubes. 

In Table |n| we give the obtained values of structural 
and mechanical properties for a set of nanotubes obtained 
from our calculations. For comparison we also give re- 
sults for the (6,6) C and BN nanotubes calculated with 
PW pseudopotential DFT calculations. As can be seen, 
the agreement between the TB and first-principles cal- 
culations in both structural and mechanical properties is 
rather good. Notice also that the values of Y s (and in 
fact those of Y also) for C, BN and BC3 nanotubes of 
similar diameters are approximately in the same ratio as 
predicted from Eq. (0). 
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FIG. 2. Young's modulus as a function of the tube diam- 
eter for C, BN, BC 3 , BC 2 N (structure II only) and C3N4, 
as calculated from the tight-binding simulations. Results ob- 
tained for (n,n) nanotubes (filled symbols), (n,0) nanotubes 
(empty symbols) and also for C (10,5) (+) and (10,7) (x) are 
shown. 




FIG. 3. Relaxed structure for the C3N4 (8,0) nanotube, as 
obtained from the TB calculations. The dark shaded atoms 
are Nitrogen atoms, while the lighter ones are Carbon atoms. 

In Fig. U the values Y s have been plotted as a function 
of the tube diameter for the different types of tubes con- 
sidered in this work. The first observation that can be 
extracted from Table || and Fig. || is the fact that car- 
bon nanotubes are predicted to have the highest Young's 
modulus of all the different types of tubes considered. BN 
and BC3 tubes have very similar values of Y s , though the 
latter have slightly larger values. The BC2N nanotubes 
are predicted to be slightly stiffer than the BN and BC3 
tubes, while C3N4 nanotubes lie well below the rest in 
stiffness. The value of Y s we obtain for the wider C nan- 
otubes, 0.43 TPa nm, corresponds to a Young's modulus 
of 1.26 TPa in the conventional definition of Eq. (|l|), if 
we take SR = 0.34 nm, i.e. the inter-layer spacing in 
graphite. This value is in very good agreement with the 
experimental value recently obtained by Krishnan and 
coworkers @ for SWNT's (1.25 TPa). It is also in rather 
good agreement with the value of 1.28 TPa reported by 



5 



Wong et al. |24| , though this later value was obtained for 
MWNT's. However, it is expected that the Young's mod- 
ulus be mostly determined by the intra-wall C-C bonds, 
and it is therefore not surprising that the values look so 
similar for both MW and SWNT's. For composite nan- 
otubes, the only experimental data on mechanical prop- 
erties currently available to our knowledge are the re- 
sults of Chopra and Zettl |23| , who have measured Y for 
BN MWNT's. They quote a value of 1.22 TPa, which 
is somewhat larger than the result we obtain for these 
tubes, but nevertheless the agreement is close. Tough it 
may seem that the choice 5R = 0.34 nm is somewhat 
arbitrary, it should be pointed out that the experimen- 
tal results are not free of this arbitrariness either, given 
that to obtain a value of the Young's modulus according 
to Eq. (|l|), it is necessary to interpret the experimen- 
tal observations on the basis of some mechanical model, 
usually a hollow cylinder with a certain wall thickness. 
Clearly, in the case of SWNT's the question of how to 
choose SR — 0.34 nm applies to experiments as well as 
to theoretical calculations. 





FIG. 4. Relaxed structure for the C3N4 (5,5) nanotube. 

The reason why C3N4 nanotubes are predicted to be 
so much softer than all other types of nanotubes consid- 
ered in this work is the fact that for a given amount of 
tube surface, tubes of this composition have a smaller 
density of chemical bonds. Indeed, these tubes present 
a more hollow structure when compared to the other, 
perfectly hexagonal nanotubes (see Figs. || and It 
is also interesting to note that, while C, BN, BC3 and 
BC2N nanotubes do not show noticeable differences for 
the structural or mechanical properties as the chiral an- 
gle is varied, sizeable differences are observed in the case 
of C3N4 nanotubes. It can be seen in Table || that the 
(n,n) C3N4 tubes have a higher Young's modulus than 
the (n,0), though both seem to converge towards the 
same number in the limit of large diameters. This dif- 
ference is also reflected in the structure, as can be seen 
in Figs. H and [|. Notice how the the hexagons are at an 
angle with the surface of the tube which is different in 
(n,n) and (n,0) tubes. 

As for comparison with other theoretical predictions, 
the results quoted by Lu 1 37 ] are somewhat smaller than 
ours. The results in ref. [37| are 0.97 TPa for all tubes. 
This difference is most likely due to the different ap- 
proaches (empirical potentials and TB) used in that work 
and ours. We also observe in our results for Y s a slight de- 
pendence on the tube diameter. As the diameter becomes 
larger, Y s approaches a plateau value which corresponds 
to the value calculated for the flat graphene-like sheet of 
each nanotube composition. Interestingly, the approach 
to the limit value is from below, as can be expected if 
one considers that bending a flat graphene sheet weakens 
the bonds. Given that it is the strength of the chemical 
bonds which determines the actual value of the Young's 
modulus, it is natural that small-diameter (high curva- 
ture) tubes have smaller Young's moduli, and in the limit 
of large diameters, the mechanical properties essentially 
correspond to those of the flat graphene sheet. In con- 
trast, the results of Lu (3^] are largely insensitive to the 
tube diameter. This is due to the fact that pair-potential 
models, such as the one used by Lu, do not reflect the 
changing nature of the chemical bonding as the curvature 
is changed. To reproduce this effect, a model sensitive to 
the changing environment (i.e. a many-body model) is 
required. 

To summarize, we have used a non-orthogonal TB 
model parametrised for C, B and N based systems to 
perform a systematic study of the energetic, structural 
and mechanical properties of single-wall nanotubes of 
different chemical composition. We have checked the 
accuracy of our predictions against some first-principles 
calculations, and the agreement obtained is good. Fur- 
thermore, we obtain good agreement with the available 
experimental data. We obtain strain energy vs. diame- 
ter curves which obey very closely the expected D~ 2 be- 
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haviour. Our results show that carbon nanotubes are ex- 
pected to be stiffer than any of the composite nanotubes 
considered, having a Young's modulus of approximately 
1.3 TPa, which corresponds to that of a flat graphene 
sheet within the same theoretical model. The C3N4 nan- 
otubes, which present a more hollow structure than the 
other tubes, are predicted to have a Young's modulus 
nearly half that of the carbon nanotubes. 
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